function [ A_ineq,b_ineq] = constraint_mat_theta_2(combos_all_unique,J)

% monotonicity

[ A_monot ] = A_monot_fun_2(combos_all_unique,J);

% [A_exch] = A_exch_fun(m_regr,J);

A_ineq = sparse([A_monot]);

b_ineq = [-ones(size(A_ineq,1),1)*0.0];

end

